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ON THE MINIMAL TRAVEL TIME NEEDED TO COLLECT n 
ITEMS ON A CIRCLE 

By Nelly Litvak and Willem R. van Zwet 

University of Twente and University of Leiden 

Consider n items located randomly on a circle of length 1. The 
locations of the items are assumed to be independent and uniformly 
distributed on [0, 1). A picker starts at point and has to collect all n 
items by moving along the circle at unit speed in either direction. In 
this paper we study the minimal travel time of the picker. We obtain 
upper bounds and analyze the exact travel time distribution. Fur- 
ther, we derive closed-form limiting results when n tends to infinity. 
We determine the behavior of the limiting distribution in a positive 
neighborhood of zero. The limiting random variable is closely related 
to exponential functionals associated with a Poisson process. These 
functionals occur in many areas and have been intensively studied in 
recent literature. 

1. Introduction. This paper is devoted to the properties of the optimal 
route of the picker who has to collect n items independently and uniformly 
distributed on a circle. By optimal route we mean the route providing the 
minimal travel time (see Figure 1). The problem has applications in perfor- 
mance analysis of carousel systems. A carousel is an automated storage and 
retrieval system which is widely used in modern warehouses. The system 
consists of a large number of shelves or drawers rotating in a closed loop in 
either direction. Orders are represented by a list of items. The list specifies 
the type and retrieval quantity of each item. The picker has a fixed posi- 
tion in front of the carousel, which rotates the required items to the picker. 
In this paper we study the minimal travel (rotation) time of the carousel 
while picking one order of n items, the locations of which are assumed to be 
independent and uniformly distributed on the carousel. 

Let Uq = be the picker's starting point and, for i = 1, 2, . . . ,n, let the 
random variable Ui denote the position of the ith item. The random vari- 
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D n +l,n * Dl.n 



Fig. 1. Minimal travel time on a circle. 

ables Ui, U2, ■ ■ ■ , U n are independent and uniformly distributed on [0,1). 
Set U n +i = 1. Let 

0=Uo:n <Ul; n < ■ ■ ■ <U n -.n < U n +\ : n = 1 

denote the ordered C/o, C/i, . . . , J7 n +i- Then the picker's starting point and 
the positions of the n items partition the circle into n + 1 uniform spacings 

Di,n = Ui :n -Ui-i :n , l<i<U + l. 

Let Xi,X2, . . . be independent exponential random variables with mean 1 
and write 

i 

S = 0, Si = J2Xj, »>i- 

i=i 

It is well known that [cf. Pyke (1965)] 

(1.1) (-Di, n , D2 >n , ■ ■ - ,D n +l,n) = (Xi/ S n+ i,X2/S n+ i, . . . , X n+ i/S n+ i), 

that is, the spacings are distributed as normalized exponentials. According 
to Pyke (1965), this construction is useful "to show that an ordering of uni- 
form spacings may be considered as an ordering of the exponential random 
variables." 

Now, let T n be the minimal travel time. We explore T n in terms of the 
uniform (n + l)-spacings D\^ n ^2,n^ ■ ■ ■ ,D n +i jn . For n = 1, the problem is 
trivial. The picker just chooses the shorter distance from the starting point 
to the item, and thus the travel time T\ is distributed as (1/2)1)^1 (a nor- 
malized minimum of two exponentials). For n = 2, one can easily verify that 
the optimal route is guaranteed by the nearest item heuristic where the next 
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item to be picked is always the nearest one. The travel time distribution for 
this heuristic was obtained by Litvak and Adan (2001). It follows from their 
result that T2 is distributed as (\/2)D\ t 2 + (3/4)I?2,2- For n > 3, the problem 
becomes much more difficult. 

A crucial and simple observation made by many authors [see, e.g., Bartholdi and Platzman 
(1986)] is that the optimal route admits at most one turn. Obviously, it is 
never optimal to cover the same segment of the circle more than twice. Thus, 
in general, T n can be expressed as 



(1.2) 



T n = 1 - max< max {Dj n - Uj-i ■ „}, 



max {D n+2 -j,n - (1 - U n+2 -j:n)} >■ 

i<j<n J 



This formula is easy to understand by means of Figure 1. Clearly, for j = 1, 
2, . . . , n, the term Dj tU — C/,_i : n is the gain in travel time (compared to one 
full rotation) obtained by skipping the spacing Dj^ n and going back instead 
(ending in a clockwise direction). The same can be said about D n +2-j,n ~ 
(1 — U n+ 2-j :n ), but here the picker ends in a counterclockwise direction. 
Under the optimal strategy, the picker chooses the largest possible gain. 

Let Tn" 1 ^ be the travel time under so-called m-step strategies: the picker 
chooses the shortest route among 2(m + 1) candidate routes that change 
direction at most once (as does the optimal route) and only do so after 
collecting no more than m items. It was proved by Litvak and Adan (2002) 
that for 2m < n, 

T^ 1 = 1 — max< max {Dj n — Uj-i ■ „,}, 
[l<j<m+l J ' J 

max {D n+ 2-j >n ~ (1 - U n+ 2-j:n)} \ 
l<j<m+l J 

/-, n\ = 1 — t, m ax| max {Xj — 

(1-3) Sn+i U<i<m+i J J 

max {X n+ 2-j — (S n+ i — S n+ 2-j)} } 

l<j<m+l J 

{m+l ^ m+l ^ "| 

23 — 1*^'"' 51 2i — i^ n+2 -3> n j ■ 



= 1 — max< 



Formula (1.3) follows from the following curious property of exponential 
random variables obtained by Litvak (2001). 

LEMMA 1.1. For any m = 0,1,... and 0<q<l, 

m+l 

(1.4) max {Xj - {q^ - l)S^} ± (q~ l - 1) £ q\l - q^Xj. 

!<]<m+l 
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The proof also implies that for any m = 0, 1, . . . , n, 

m+l 

(1.5) max \D hn - {q~ l - 1)E7,_ 1:B } = (g" 1 - 1) £ g*(l - <f) l ^>- 

l<7<m+l — ; 

J=l 

If 2m < n, then the two internal maxima in the third expression of (1.3) 
are independent and (1.4) can be used to rewrite each of them separately. 
Moreover, the same argument applies for the normalized exponentials yield- 
ing (1.3). If 2m > n, then the two internal maxima become dependent and 
the argument fails. 

In fact, the optimal strategy is the m-step strategy with m = n — 1. In- 
tuitively, it is clear, however, that with high probability the optimal route 
has only a few steps before a turn. That is, the m-step strategy often pre- 
scribes the optimal picking sequence even when m is relatively small. It was 
shown by Litvak and Adan (2002) that already for m = 2, the m-step strat- 
egy is quite close to optimal and, on average, outperforms the nearest item 
heuristic. 

Let Kn"^ and K n denote a number of items collected before a turn un- 
der the m-step strategy and the optimal strategy, respectively. If there is no 
turn, these numbers are set equal to zero. It was proved by Litvak and Adan 
(2002) that: 

(i) Tr™^ and are independent random variables; (ii) for any k = 0, 1, . . . ,m 
and 2m < n, 

P(i^M = k) = F[ [arg max {D in -U i - 1 . n \ = k + l 

(1.6) = pf[arg max {X 7 - — Sj-i} = k + l 

\l I<i<™+1 
1 



<2k—m ' 

(hi) for any k = 0, 1, . . . , n — 2, 

(1.7) ¥(K n >k)< l/2 k . 

The last estimate is helpful in the analysis of the limiting properties of the 
optimal route. For example, it was proved by Litvak and Adan (2002) that 
for any fixed k = 0, 1, . . . , 

(1.8) lim F(K n = k) = l/2 k+1 . 

n— >oo 

Indeed, observe that for k = 0,1, . . . ,m, 

P(i^M = k) - F(K n > m) < F(K n = k)< F{K^ = k). 

Now, let m and n go to infinity in such a way that the inequality 2m < n is 
always satisfied. Then (1.8) follows readily from (1.6) and (1.7). 

In this paper we first derive simple upper bounds for the minimal travel 
time. Then we analyse the distribution of T n . Further, we obtain the limiting 
behavior of T n when n tends to infinity. 
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2. Upper bounds. Let T n be the minimal travel time needed to collect n 
items independently and uniformly distributed on a circle of length 1. The 
following lemma gives an upper bound that holds for any realization of the 
random items' locations. 

Lemma 2.1. For any n > 1, the travel time T n never exceeds 1 — a n +i; 
where 

1 

an+i = —, , n = 2m: 

K-r-L 2«i+i _|_ 2 m — 2 

1 

ctn+i = n 1 n = 2m+l. 

' 2 • 2 m + — 2 

TMs upper bound is tight. 

Proof. Assume that n = 2m + 1. For n = 2m the proof is similar. The 
positions of the items plus the picker's starting point partition the circle 
into n + 1 spacings with lengths d\ , d 2 , . . . , . Note that for any collec- 
tion di,d2, ■ ■ ■ , d n +i > there exists a number j = 1, 2, . . . ,m + 1 such that 
either (i) dj > 2 j ~ 1 a n+ i, di < 2 l ~ 1 a n+1 , I = 1, 2, . . . , j - 1, or (ii) d n+2 -j > 
2 jf ~ 1 a n+ i, d n+ 2-i < 2 l ~ 1 a n+ i, I = 1, 2, . . . , j — 1. This follows since 

m+l 

2 ^ 2 J '- 1 a n+ i = di + d 2 + ■ ■ ■ + d n+1 = 1. 
j'=i 

Without loss of generality assume (i). Then the route that skips the spacing 
dj and goes back instead has length 

1 — dj + d\ + d 2 H h dj-i < 1 — a n+ i , 

and its length must be greater or equal than T n . This proves the upper 
bound. 

To show the tightness we just put dj = d n+ 2-j = 2 J_1 a n +i , j = 1, 2, . . . , m + 1. 
In this case the travel time under the optimal strategy equals 1 — a n +i- d 

Let us now consider the following approximation of T n in (1.2), 

Z? = 1 — max! max {X, — 

5" n+ i ll<i<m+l L J J 

.max , - (^n+l - S n+ 2-j)} (, 

l<j<m'+l J 

where m = m' = (n — l)/2 if n is odd and m = m' + 1 = n/2 if n is even. In 
both cases we have m + m' = n — 1 so that the -Xj's from the first internal 
maximum are not involved in the second internal maximum. That is, the 
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two internal maxima are independent, and we can apply Lemma 1.1 to each 
of these separately to arrive at 

{m+1 i m'+l °\ 

Clearly, T° gives a tight stochastic upper bound for T n . In fact, T° and T n 
differ with probability of order 2~ n / 2 according to (1.7). It was shown by 
Litvak and Adan (2002) that T° is stochastically larger than the weighted 
sum 



n+1 



T n = E( 1_a i) jD J>- 

Straightforward estimation of the expected difference between T* and T° 
yields 

E(7;°-T n *)<0.09E(ZV) 



n + 1 



Thus, 



1 ™±} s 0.09 



(2.2) E(T n ) < E(T n °) < — 2(1 - «i) + n 

In Table 1 (see Section 4), we compare the mean travel time E (T n ) obtained 
by simulation with upper estimate (2.2) and approximation (4.8), which 
follows from the limiting results in Section 4. The results prove that the 
bound (2.2) is quite sharp. For larger n, however, (4.8) gives a slightly better 
approximation. 



3. The minimal travel time distribution. In this section we produce an 
explicit expression for IP(T n > 1 — t). First, note that it is never optimal to 
turn after covering half of a circle. Now, consider the events 

A n>k (u,v) = [Uk:n =U < 1/2 < 1 - V = U k+1:n ], 

0<«,u<l/2,Jb = 0,l,...,n. 

For k = 2, 3, . . . , n — 2, the joint distribution of Ui :n , . . . , £/fe-i :n; 1 — Uk+2-.m ■ 
1 — U n - n given A n ^{u,v) is that of 



vHl:k-l, • • • i u Uk-l:k-l, v Vn-k-l:n-k-li ■ ■ ■ , v Vi; n -k-l> 
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where U and V are independent vectors of uniform order statistics. As the 
event [T n > 1 — t] implies 1 — v — u — uAv<t, we have for k = 2, 3, . . . , n — 2, 

¥(T n >l-t\A n>k (u,v)) 

(3.1) 



P max {(C/ j:fc _i - Uj-uk-i) - Uj-uk-i} <t/u, 
i<j<fe 



max {(V,- :n _ fc _i - K-_i :n _ fe _i) - K_i :n _ fe _i} < t/v 

l<]<n—k 
x v— u— uAv<t] 

[1— v— u— uAv<t] ■ 

Here u A v = min{ii, v} denotes the smaller of u and v and 



(3.2) P m {t)=¥( max {D j>m - U^ Um } < t , m = 1,2, . . . ; t > 0. 

\l<j<m+l / 

One readily verifies that the final expression in (3.1) continues to hold for 
k = 1 and k = n— 1 , provided we define 

(3.3) P (i) = l [t>1] , t>0. 
For k = and k = n, we find 

P(T n > 1 - t\A nfi (0,v)) = P n _i(i/^)l[i-„< t] , 
P(T n > 1 - t|A,n(n, 0)) = P n _l(V!i)l[l- u < t] . 

It follows that 



c (T n > 1 - t) 

r-1/2 /-1/2"-- 1 



(3.4) 



o JO 



k=l ^ ' 



n— fe— 1 



X f5fc-l(t/'")- P n-fc-l(*/' i; ) 1 [l-^-«-«A«<t] ^cfo 
/•1/2 _ 

+ 2 • l[ t>1 / 2 ] / l P n -i(t/u) du. 

Ju=l—t 



Formula (1.5) and Theorem 2 of Ali and Obaidullah (1982) imply an ex- 
pression for P m (t). Writing 

c j = (2P-l)- 1 , J = 1,2,..., 

and x+ = max{x, 0} for the positive part of a number x, we find that for m = 
1, 

2,..., 

(m+l \ rn+1 m+1 

£ Cj ^, m < * = £ {(* - Cj )+r n fa - c i) _1 - 
j=l / j=l (=1, 
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The last expression is also valid for m = 0. Of course, for t > 1, the terms in 
(3.5) sum to 1. 

Alternatively, one can determine P m (t), recursively. Conditioning on Uxirm 
we find for m = 2, 3, . . . , 



P[ max {D im -Ui-i :m }<t Ui :m = u) 

\l<%<m+l J 



= P[(l-u) max {A,m-1 ~ A-l: m-l} ~ « < < lru<tl 
V l<i<m / — 

This yields the recursive equation 

(3.6) P m (t) = jf* m(l - n)™" 1 ^-! (^) d«, 

which is valid for m = 1,2, 

We can now find the distribution of T n by substituting (3.3) and either 
(3.5) or (3.6) in (3.4) and integrating. One obtains, for example, for t > 0, 

Pi(*) = i(3t-l) + -§(t-l) +J 

p 2 (*) = H(w - 1) + } 2 - im - 1) + } 2 + |{(t - 1) + } 2 

and for < i < 1, 

P(Ti>l-t) = (2t-l)+, 

P(T 2 > 1 - 1) = \{{At - 1) + } 2 - 2{(2i - 1) + } 2 , 

P(T 3 > 1 - 1) = i{(6t - 1) + } 3 - §{(4t - 1) + } 3 

-|{(4t-l) + } 2 + ^{(2t-l) + } 3 . 

Although the general structure of these functions is fairly easy to under- 
stand, it seems quite useless to provide explicit expressions for P(T n > 1 — t) 
for much larger values of n. Instead, we study their asymptotic behavior in 
Section 4. 



4. Limiting results. In this section we shall obtain the limiting distri- 
bution of (n + 1)(1 — T n ). First of all, let us consider the limiting behavior 
of 



P m (t/(m+l)) 



( m+1 1 \ 

P^( ? n+l)^-- T D J , m <tJ. 
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Theorem 4.1. Let X\,X2, ■ ■ ■ be independent exponential random vari- 
ables with mean 1. Then 

m+l -. oo -. 

(4.1) (m + 1) £ -- -D jtm 4 £ — -A- 

and the limiting distribution is given by 
P(t) = lim P m (t/(m + l)) 

m — >oo 

1 ' J =l-^(-ir 1 2^exp{-(2^-l)i}n^r L T- 

3=1 1=1 

The distribution function P satisfies the integral equation 

t-2t 



(4.3) e~ t P(t) = J e- u P(u) 



du. 



Proof. The argument essentially repeats the proof of Theorem 4 of 
Litvak and Adan (2001). Define 

m+l -. oo -. 

(4.4) j m = J2 ajri^' J = E 2TTi x i- 

i=i i=i 

By the monotone convergence theorem, E(</) — \i'ni rn ^. Q E( J m ) < oo. In par- 
ticular, it implies P(J < oo) = 1. 
Now, using (1.1), we write 

t i_ i\ V" 1 n d (m+ l)J m 
(m + 1) 2^ - — -D j>m = — . 

j=1 U — 1 Jm+l 

By definition, the sequence {J m } converges a.s. to J. The strong law of 
large numbers implies that the sequence {(m + l)/S m+ i} converges a.s. to 
1. Thus, {(m + 1) J m /S m+ i} converges a.s. to J which immediately gives 
(4.1). 

The distribution P of J can be obtained via inversion of its Laplace- 
Stieltjes transform 

SSL. 2? — 1 
<p(s) = E(exp(-aJ)) = J] 2 /_ 

3=1 

One can expand <p(s) in rational fractions of s and obtain 

°° (-1)3-103 Ci 1 
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Here, in order to write the formula for the residues of tp(s), one can apply 
well-known expressions from so-called g-calculus [see, e.g., Gasper and Rahman 
(1990)], but in our case it is not difficult to verify this formula directly. In- 
version of (4.5) yields (4.2). 

Finally, we use (3.6) and the dominated convergence theorem to obtain 



m— >oo 



P(t) = Jim_P m (t/(m + l)) 

V l-u 



*/("•+!) . .„ ! ft/(m + l) + u\ 
lim / mil - u) m - l P m _i ( M . ; ) du 



m—*oo 







t e- u P{t + u)du 
Jo 

e l J e~ u P(u)du, 



which proves (4.3). □ 

Obviously, we also have convergence of moments. For the kth moment of 
P, (4.2) yields 

poo °° oj I 1 

E ( j') = / o t » (B . w = Hg(-iy-.^_ I5r n ? _. 

Alternatively, one can directly use (4.4) to find a simple expression for 
cumulants n u , v > 1, of P. It is immediate that 

oo 

e(j)=k 1 =^(2^-i)- 1 , 

oo 

Var(J) = K 2 = ^(2^-l)- 2 . 
j=l 

Furthermore, 

oo oo OO / - 1^ 

log(E exp(itJ)) =-J2 log(l - (2 J - l)" 1 *^ = ^ E 



where i is the imaginary unit. Since log(Eexp(it J)) = YlvLi K jy (it) !y (^!) , it 
follows that 

oo 

^ = ^-l)!^'-l)-, i/>l. 
i=i 

The distribution function P on [0, oo) has the remarkable property that 
it is infinitely often differentiable and that all of its derivatives P^ vanish 
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at the origin. This is most easily seen by differentiating (4.3), but one may 
also use (4.2) to show analytically that p( k \0) = for all k = 1, 2, ... . It 
follows that P is not analytic at the origin. The series (4.2) diverges for all 
t < and, hence, P cannot be represented by its Taylor series around t = 0. 

Now repeating the argument from the proof of Theorem 4.1, one can show 
that 

{OO -. OO -. ^ 

where X%, X 2 , ■ ■ ■ , X[, X 2 , ■ . ■ are independent exponentials with mean 1. 
Since the two sums in the maximum are independent and (1.7) ensures 
that 

P(T n /T°) <2-( n ~ 2 )/ 2 , 
we have proved the following statement. 

Theorem 4.2. Let X\,X 2 , . . . ,X[,X 2 , . . . be independent exponential 
random variables with mean 1. Then 

{00 1 00 1 ~> 

i=i 3=1 J 

and the limiting distribution is 

(4.7) lim F(T n > 1 - t/(n + 1)) = [P(t)\ 2 , 

n — >oo 

where P(t) is given by (4.2). 

Again we have moment convergence and for the kih moment we find 
lim E[(n+l)(l-T n )] fc 

n— >oo 

00 2 j j 1 

1 (23 - 1)* l\ 2l - 1 



00 00 oi+j 3 1 * — 1 1 

- 2k\ y Y(-iy+j——L — _ fy _J_ rr _ !_. 

7=1 i=l v ' (=1 r=l 



2 i+ ' J ' 

j=l i=l 

An equivalent expression for the expectation can be obtained as 
lim E[(n+l)(l-T n )] 



(l-[P(t)] 2 )dt 
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oo -, oo oo ni+j 3 i ' 1 

= 2 y^— -yv(-i)^— — tt —r- — tt — - — 

j=l i=l j=l (=1 r=l 

« 2.1578. 

For large n we therefore have the estimate 

2.1578 



(4.8) E(T n ) ~ 1 



n + 1 



In Table 1 we compare the mean travel time obtained by simulation with 
upper estimate (2.2) (see Section 2) and approximation (4.8). We see that 



Table 1 

Estimation of the mean travel time under the optimal strategy 



n 


3 


5 


10 


15 


20 


30 


E(T n ) 


0.5262 


0.6591 


0.8052 


0.8653 


0.8972 


0.9304 


E[(ra + l)(l-T n )] 


1.8952 


2.0454 


2.1423 


2.1548 


2.1592 


2.1572 


Upper estimate (2.2) 


0.5433 


0.6670 


0.8068 


0.8658 


0.8976 


0.9306 


(n + l)[l-upper estimate(2.2)] 


1.8268 


1.9980 


2.1252 


2.1472 


2.1504 


2.1514 


Approximation (4.8) 


0.4605 


0.6404 


0.8038 


0.8651 


0.8972 


0.9304 


(n + l)[l-approximation (4.8)] 


2.1578 


2.1578 


2.1578 


2.1578 


2.1578 


2.1578 



both approximations are quite sharp, but (4.8) performs somewhat better. 
It is no surprise that both (2.2) and (4.8) are close to E(T n ) for large n 
since all three quantities converge to 1 as n — > oo. What is encouraging is 
that, already for n = 30, both approximations of (n + 1)(1 — E(T n )) are very 
good. That (4.8) yields a better approximation of E(T n ) than (2.2) is to be 
expected since it is asymptotically correct up to and including order n _1 , 
whereas (2.2) has a slight asymptotic error of about + 0.006/(n + 1). After 
all, (2.2) was derived as an upper bound. 

5. Asymptotic behavior in the neighborhood of zero. In this section we 
study the behavior of P(t) as t — > +0. So far we have found only that P 
has vanishing derivatives at the origin and can not be expanded in a Taylor 
expansion around t = 0. We shall, therefore, have to attack this problem in 
a different manner. 

Let Xi,X2,.-. be independent exponential random variables with mean 
1, let 



Cj = {2 j - I)" 1 , J = 1)2, ... , 
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and define 

oo 

J = H c i X r 
3=1 

We want to determine the behavior of 

P(t) = P(J<t) 

for small positive values of t. In principle this problem is solved in Theo- 
rem 3.2 of Davis and Resnick (1991), but we need to do a substantial amount 
of analysis to make their result explicit, even in our relatively simple case. 

In our case, the distribution function F(x) = P(X% < x) = 1 — exp{— x} 
and the density f{x) = exp{— x} are regularly varying at with index a = 1 
and a — 1 = 0, respectively. The c,-'s are positive and nonincr easing, their 
sum converges and for every 9 £ (0, 1), 

oo oo 

9 n J2{c 2 3 /c 2 n }i [j>e - n] =e n Y,{(2 n - i)/(2 j - 1)} 2 W"] -o 

3=1 3=1 

as n — > oo. The density / satisfies 

;>oo 

/ e- 2Xx f 2 (x)dx = l/{2(l + A)} forA>0. 
J o 

Hence, we have verified the assumptions of Theorem 3.2 of Davis and Resnick 
(1991) in our case. The theorem states that 



(5.1) P(m\) ~ exp{Xmx}(fj(X)/(\SxV2TT ) asA^cx). 

Here 

OO CO ^ 

mA= gi^ = g273ITA' 

OO ^ OO nj ^ 

= n i+j^ = n 2i3i+A 

3=1 J 3=1 



and 



S 



c 2 

2 — 3 = \ "* x 

^(1 + Ac,)2 2^ {2J _ 1 + X y 



A 



3 = 

We obviously have to study the behavior of these quantities as A — > oo and, 
hence, m\ — > 0. It is easier to deal with integrals than sums. For k = 1, 2, . . . 
and A — ► oo , we have 

0< / (2 X - 1 + X)~ k dx - V(2 J '-l + A)~ fc 
Jo 3=1 
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< £(2^-1 + A)" fc 1 + A)- 
3=0 i=i 

= A~ fc 



and, hence, 



j 

For k = 2, this yields 



°° POO 

- 1 + A)~ fe = / (2 X - 1 + A)" fc dx + 

7 = 1 ^° 



0(A 



roo 

(52) 5^ = aog2)- 1 y (y + A)- 2 (y+l)- 1 ^ + 0(A' 

= (logA)/(A^log2) + 0(A- 2 ), 



as A — > oo. If we apply the same approach to Am^ and logyj(A), however, 
then the error caused by approximating these sums by integrals is of the 
order 0(1) and O(logA), respectively, which yields a multiplicative error 
factor (1 + 0(a\ b )) in (5.1) for some positive a and b. Of course this is not 
good enough so we shall have to expand the series representing \m\ and 
logy j (A) directly with remainder o(l) in both cases. 
Let A; be a natural number and 9 G [0, 1) be such that 

A = 2 fe+e , 

and thus 

k = (log A)/ log 2 — 9 = [(log A)/ log 2J , 

9 = (log A)/ log 2 - k = frac((log A) / log 2) . 

Here [x\ and frac(x) are the integer and the fractional part of x, respectively. 
In order for A — > oo, it is necessary and sufficient that k — > oo, while 9 may 
vary arbitrarily in [0, 1) with k. Using (5.2) we find 

oo ^ 

Am A = 



i=1 »-l + A 
E 23-1 + 2 k+e 

3=1 
oo 2 fc + 6 ' 

: E27T2^ + 0(A_ll ° gA) 

k oo 

E^TT+.E 2?=^ + OCA" 1 log A) 
j=i j=fc+i 
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k— 1 -. oo -. 

3=0 3=1 
k Qp oo ^ 

= E + g 2F^TT + 0(A_1 log A) 

fc 9 1_e 00 2 d 
= logA_ 2 + y _j_^_g + Q(A -l 1 A) , 

Hence, 

(5.3) Am A = (log A)/ log 2 + A(0) + 0(A~ 1 log A) , 



with 



oo 2 1_e 00 2 6 



A (°) ~ ~ E 2 , + 21=5 + E 27T¥ " °- 

3=1 3=1 

Notice that the term ^4(#) of order 1 is not constant but depends on 9 S [0, 1). 
The expansion (5.3) is obviously uniform in 9. 
Similarly, by (5.3), 

oo 

log^(A) = X>g{(2 j - 1)/(2 J - 1 + A)} 

3=1 

oo oo 

= 5>g{2V(2>' + A)} + £log(l - 2~ J ) 

3=1 3=1 

oo 

+ 5>g{l + 1/(2^-1 + A)} 

3=1 

oo oo 

= Y, log{2V(2 J + 2 k+e )} + £ log(l - 2-*) 

3=1 3=1 

+ 0(A~ 1 logA). 

Furthermore, 

oo 

5>g{27(2^ + 2 fc+9 )} 

3=1 

fc oo 

= £ log{2 1 - 9 /(2^ + 2 1 "*)} + £ log{27(2* + 2 e )} 

3=1 3=1 
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= £(1-0) log 2- (l/2)fc(fc + l)log2 

k oo 

- e iog(i + ^~ e ~ 3 ) - E Mi + 29 " J )- 

3=1 J=l 

Substituting A; = (log A) / log 2 — 6 and using 

oo oo 

E log(l + 2 1 ^)< E 2 1 -^ = 0(A" 1 ), 

j=fc+l j'=fc+l 

we finally find 

(5.4) log^(A) = + ^ + 5(0) + OCA" 1 log A), 

where 

oo 

<r « B W = E lQ g{(l " 2" J )/[(1 + 2^')(1 + 2 1 -^')]} 

(o-oj i=1 

-(1/2)0(1-0) log 2. 

Again the term B{9) of order 1 depends on 9 and the expansion is uniform 
in 0G [0,1). 

Substituting (5.3), (5.4) and (5.2) in (5.1), we obtain, for A — > oo, 



(5.6) 



1 1 

+ I 77 + 



logA + A(0) + 5(0)1. 



,2 log 2/ v ' v ') 

It remains to find approximations of log A and (log A) 2 as functions of 
t = m x = (log A)/(A log 2) + A{9)/\ + 0(A^ 2 log A). 

We find 

log(l/i) = log A - log log A + log log 2- A(9) (log 2) / log A 
+ 0((logA)- 2 ), 
loglog(l/i) = log log A - (log log A) / log A + (log log 2)/ log A 
+ 0((loglogA) 2 /(logA) 2 ) 

and, hence, 

log A = log(l/t) + loglog(l/t) - log log 2 + (loglog(l/t))/log(l/t) 
- (loglog2)/log(l/t) + A(0)(log2)/log(l/t) 
+ 0((loglog(l/t)) 2 /(log(l/t)) 2 ), 
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(log A) 2 = pog(l/i) + log log( 1/t ) - log log 2] 2 + 2 log log (1/t) 
- 2 log log 2 + 2A(9) log 2 
+ 0((loglog(l/t)) 2 /log(l/i)). 

Together with (5.5) and (5.6), this yields that, for t — > 0, 

P(t)~C(0)exp{-(21og2)" 1 [log(l/t)+loglog(l/t)-loglog2] 2 } 

x r (l/2+l/log2) 

with 

(5.7) c( ^ 2 -.a-.)/^n (1+2 4)" (1 2 ;V.- Jr 

The factor C(6>) depends on 6* = frac((logA)/log2). 
It remains to express 9 in terms of t. Define 

i>(t) = (log2)- 1 [log(lA) +loglog(l/t) - log log 2]. 

We have k + 6 = (log A)/ log 2 = ip(t) + o(l), and as C is positive and bounded, 
the derivative of C is positive and bounded and C(9) = C(1 — 6). This im- 
plies that C(frac{V>(*)}) = C(0)(1 + o(l)). It follows that, as i -> +0, 

(5.8) P(t)~C(frac{^(t)})exp(-^[?/;(t)] 2 lt-( 1 / 2+1 / log2 ), 



with C defined in (5.7). This is an exact asymptotic expression for P(t) as 

The dependence on hac(^(t)) in (5.8) is a most unusual feature. In fact, 
preliminary numerical calculations make one wonder whether there is any 
dependence at all, since one finds that C(9) equals a constant (~ 0.01013) 
throughout the interval < 9 < 1 to any reasonable degree of accuracy. Thus, 
in order to properly understand the asymptotic expression (5.8), we have to 
analyze C(9) in more detail. Proposition 5.1 states that C{9) does indeed 
depend on 9, but in a very peculiar way. In fact, for any real 9, 



C{9) 

where 



2V827T 



{h{0)y l ~ 0.01013(^3 (A)) -1 , 



K Q\ ^ = 1 + 2 E exp{-2fc 2 7r 2 / log 2} cos{2fc7r(l/2 - 8)} 

{ ' k=l 

= 3 (tt(1/2 - 9), exp{-2vr 2 / log 2}). 
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Here $3 is a theta function 

00 

& 3 (z, q) = 1 + 2 ^ q k ' 2 cos(2kz) . 
fc=l 

Note that for all 9, 

\& 3 (P)-i\ <ht 12 

is a quantity which is difficult to reveal numerically! 
Proposition 5.1. For any real 9, 

n(i+2 e - j ')(i+2 1 ^) 

(5-10) ol /g nr- 00 

=2-^)/^ 3 w^^ri( i -2- j )- i 5 

Vlog2 

where $3(0) is given by (5.9). 

Proof. We first apply Jacobi's triple product identity [see, e.g., Askey 
(1980) and Gasper and Rahman (1990)]. For any q£ (0,1), 

(5.11) J(l-xq j )(l-x" 1 q j+1 )(l-q j+1 )= £ (-l) n q®x n . 

j=0 n=—oo 

Take x = -2~ e , q = 1/2. Then (5.11) becomes 

00 00 

(5.12) 7(l + 2 e ^')(l + 2 1 - e ^')(l-2^)= 2~< n ~ 1 V 2 2~ 9n . 

j=l n=—oo 

The right-hand side of (5.12) is of the form 

00 

40) £ g(n), 

n=— 00 

where 

V ' v^g2 



g(x) = ^|Eexp{-(l/2)(log2)(x - 1/2 + #) 2 } 

V Z7T 



and 



is a normal density with mean fi = 1/2 — 9 and standard deviation c 
1 / -y/log 2. The characteristic function of g is given by 

7 (t) = exp{-t 2 /(21og2) + it(l/2 - 0)}, 
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where i is the imaginary unit. For each fixed A and for each real £, the 
Poisson summation formula [see Feller (1970)] gives 

+00 +00 
(5.13) Yl l(Z + 2k\) = j ]T <7(mr/A)exp{m(7r/A)£}. 

k=-oo n=-oo 

Put A = 7r, £ = 0. Then the right-hand side of (5.13) becomes Y^=-oo9i n ) 
and on the left-hand side we have 
00 00 

7(2/cvr) = 7(0) + exp{-2/eV/log2}exp{i(l/2-0)2fc7r} 

k=— 00 fc=-oo 

OO 

= 1 + 2^exp{-2/cV/log2}cos{2/br(l/2 _ #)} 
k=l 

Hence, (5.13) reduces to 

+00 

n=—oo 

implying that the right-hand side of (5.12) equals c(0)"&^{9). This immedi- 
ately yields (5.10). The proposition is proved. □ 

We summarize our findings in the following theorem. 

Theorem 5.2. Let X\,X2, ■ ■ ■ be independent exponential random vari- 
ables with mean 1, and let 

OO 

j=j2^-i)- 1 x ] . 



Then 



p( j < t) ~ , y- 1 - 

2 1 /82^ 3 (frac{V'(i)}) 



x exp |_^[ V ,(t)]2 j t -(i/2+i/io g 2) as t ^ +o, 
where 

i>{t) = (log2)" 1 [log(l/i) +loglog(l/t) - log log 2] 
and t?3 is defined in (5.9). 
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6. Related results. In a similar fashion we can also analyze more general 
linear combinations of i.i.d. exponential random variables than J. For any 
q G (0, 1), define 

oo 

jw = («- 1 -i)E(9" i - 1 )" ljr i- 

i=i 

Clearly, J=jW 2 \ One can show that 

m+l -, 

(m + lXg- 1 -!) E — T-{ J{<1) asm^oo, 

j=i Q 1 

where the expression on the left occurs in the right-hand side of (1.5) for 
n = m. The random variable j( q > can be written in the following way. Let 
N(t) be a standard Poisson process. Then 

(g-l _ 1)-1 J(9) = ^(*)+l(l - dt 

J 

°° /-oo 



3=1 



where 

/(<?) 



is an exponential functional associated with a Poisson process. The func- 
tional 1^ has been intensively studied in recent literature. Its density was 
obtained independently by Dumas, Guillemin and Robert (2002), Bertoin, Biane and Yor 
(2002) and Litvak and Adan (2001), for q = 1/2. Carmona, Petit and Yor 
(1997) derived a density of / °° h(N(t)) dt for a large class of functions h : N — ► 
]R+, in particular, for h(n) = q n . Bertoin, Biane and Yor (2002) found the 
fractional moments of /W. If i^ 9 '(t) is a density of JW, then i^ q \t) and all 
its derivatives equal at the point t = 0. This implies, by the way, that 
all moments of 1/1^ are finite. However, for q = l/e, it was proved by 
Bertoin and Yor (2002a) that l/I^ 1 /^ is not determined by its moments. 

The functional JW appears in a number of applications. Let T^ 1 be 
the travel time needed to collect n items independently and uniformly dis- 
tributed on a circle of length 1 operating under the nearest item heuristic 
(the picker always travels to the nearest item to be retrieved). Then it was 
shown by Litvak and Adan (2001) that (n + 1)(1 - T^ 1 ) converges in dis- 
tribution to i^ 1 / 2 ). Dumas, Guillemin and Robert (2002) showed that the 
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distribution of 1^ plays a key role in the analysis of limiting behavior of a 
Transmission Control Protocol connection. These results were extended by 
Guillemin, Robert and Zwart (2002), who found the distribution and the 
fractional moments of the exponential functional 

POO 

(6.1) i(0= e-e®dt, 



where (£(t),t > 0) is a compound Poisson process. An exponential func- 
tional (6.1) associated with a Levy process £(i) appears in mathematical fi- 
nance and many other fields. It has been studied recently by Bertoin and Yor 
(2001, 2002a, b), Bertoin, Biane and Yor (2002), Carmona, Petit and Yor 
(1997) and Yor (2001). 

Along the same lines as in Section 5, one can prove theorems similar to 
Theorem 5.2 for JW) and j( q \ In fact, it is straightforward to repeat the 
calculations for 

oo oo 

We obtain, for q G (0, 1) as t — > +0, 



/(9)<t)~ L q Vyi g(i/ g ) 



2?r 

x exp 



new) 

■i=i 



-(l/2+l/log(l/g)) 



l5^M[^)(t)]»}[^(frac{^)(t)})]-\ 

)oo 
~P(g/W<t)TT(l-^\ 

where 

^ (9) (t) = (log(l/ (7 ))- 1 [log(l/t) +loglog(l/t) - log(log(l/g))], 

oo 

#$\0) = 1 + 2^2 exp{-2fcV/k)g(l/g)} cos{2fcvr(l/2 - 0)}. 
k=l 

This agrees with the result of Bertoin and Yor (2002a) that 

logi(t) ~-i(log(l/t)) 2 ast^+0, 
where i(t) is a density of 



poo , , °° 

J ° .7=1 
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For the functional i^ 1 / 2 ) , which describes the limiting behavior of the travel 
time under the nearest item heuristic, we find 

oo 

P(/(V2) <2t)~P(J<t) Y[(l-2- j )- 1 , 

j=i 

Acknowledgments. We are grateful to Bert Zwart for drawing our atten- 
tion to relevant literature and particularly to the paper of Davis and Resnick 
(1991), which inspired the results in Section 5. Discussions with Fred Steutel 
also proved very helpful. 

REFERENCES 

Ali, M. M. and Obaidullah, M. (1982). Distribution of linear combination of exponential 

variates. Comm. Statist. Theory Methods 11 1453-1463. MR662896 
Askey, R. (1980). Ramanujan's extensions of the gamma and beta functions. Amer. Math. 

Monthly 87 346-359. MR567718 
Bartholdi, J. J. Ill and Platzman, L. K. (1986). Retrieval strategies for a carousel 

conveyor. HE Transactions 18 166-173. 
Bertoin, J. and Yor, M. (2001). On subordinators, self-similar Markov processes and 

some factorizations of the exponential variable. Electron. Comm. Probab. 6 95-106. 

MR1871698 

Bertoin, J. and Yor, M. (2002a). On the entire moments of self-similar Markov processes 
and exponential functionals of Levy processes. Ann. Fac. Sci. Toulouse Math. (6) 11 
33-45. MR1986381 

Bertoin, J. and Yor, M. (2002b). The entrance laws of self-similar Markov processes 
and exponential functionals of Levy processes. Potential Anal. 17 389-400. MR1918243 

Bertoin, J., Biane, P. and Yor, M. (2002). Poissonian exponential functionals, g-series, 
g-integrals, and the moment problem for log-normal distributions. Technical Report 
PMA-705, Univ. Paris 6, Laboratoire de Probabilites. 

Carmona, P., Petit, F. and Yor, M. (1997). On the distribution and asymptotic results 
for exponential functionals of Levy processes. In Exponential Functionals and Principal 
Values Related to Brownian Motion (M. Yor, ed.) 73-130. Iberoamericana, Madrid. 
MR1648657 

Davis, R. A. and Resnick, S. I. (1991). Extremes of moving averages of random variables 

with finite endpoint. Ann. Probab. 19 312-328. MR1085338 
Dumas, V., Guillemin, F. and Robert, Ph. (2002). A Markovian analysis of 

additive-increase multiplicative-decrease algorithms. Adv. in Appl. Probab. 34 85-111. 

MR1895332 

Feller, W. (1970). An Introduction to Probability Theory and Its Applications II. Wiley, 
London. 

Gasper, G. and Rahman, M. (1990). Basic Hypergeometric Series. Cambridge Univ. 
Press. MR1052153 

Guillemin, F., Robert, Ph. and Zwart, B. (2002). AIMD algorithms and exponential 

functionals. Technical Report 4447, INRIA. 
Litvak, N. and Adan, I. (2001). The travel time in carousel systems under the nearest 

item heuristic. J. Appl. Probab. 38 45-54. MR1816112 
Litvak, N. and Adan, I. (2002). On a class of order pick strategies in paternosters. Oper. 

Res. Lett. 30 377-386. MR1939740 



COLLECTING ITEMS ON A CIRCLE 



23 



Litvak, N. (2001). Some peculiarities of exponential random variables. J. Appl. Probab. 38 
787-792. MR1860217 

Pyke, R. (1965). Spacings. J. R. Stat. Soc. Ser. B Stat. Methodol. 27 395-449. MR216622 
Yor, M. (2001). Exponential Functionals of Brownian Motion and Related Processes. 
Springer, Berlin. MR1854494 



Faculty of Electrical Engineering 

Mathematics and Computer Science 
University of Twente 
P.O. Box 217 
7500 AE Enschede 
The Netherlands 
e-mail: n.litvak@math.utwcntc.nl 



Department of Mathematics 

University of Leiden 
P.O. Box 9512 
2300 RA Leiden 
The Netherlands 
E-MAIL: vanzwet@math.lcidenuniv.nl 



